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Abstract. The combination of dispersion measures of pulsars, distances from the model of Cordes & Lazio 1 2002 1 and 
emission measures from the WHAM survey enabled a statistical study of electron densities and filling factors of the diffuse 
ionized gas (DIG) in the Milky Way. The emission measures were corrected for absorption and contributions from beyond 
the pulsar distance. For a sample of 157 pulsars at \b\ > 5° and 60° < I < 360°, located in mainly interarm regions within 
about 3 kpc from the Sun, we find that: (1) The average volume filling factor along the line of sight f v and the mean density 
in ionized clouds n c are inversely correlated: J v (n c ) = (0.0184±0.0011)n c ~ 107±0 03 for the ranges 0.03 < n c < 2cm" 3 
and 0.8 > f v > 0.01. This relationship is very tight. The inverse correlation of f v and n c causes the well-known constancy 
of the average electron density along the line of sight. As f v (z) increases with distance from the Galactic plane \z\, the 
average size of the ionized clouds increases with \z\. (2) For \z\ < 0.9 kpc the local density in clouds n c (z) and local 
filling factor f(z) are inversely correlated because the local electron density n c (z) = f(z)n c (z) is constant. We suggest 
that f(z) reaches a maximum value of > 0.3 near \z\ — 0.9 kpc, whereas n c (z) continues to decrease to higher \z\, thus 
causing the observed flattening in the distribution of dispersion measures perpendicular to the Galactic plane above this 
height. (3) For \z\ < 0.9kpc the local distributions n c (z), f(z) and rig(z) have the same scale height which is in the range 
250 < h ^ 500 pc. (4) The average degree of ionization of the warm atomic gas J w {z) increases towards higher | z | similarly 
to J v (z). Towards \z\ = lkpc,/ v 0) = 0.24 ± 0.05 and 7 w (z) = 0.24 ± 0.02. Near \z\ = 1 kpc most of the warm, atomic 
hydrogen is ionized. 
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1. Introduction 

The interstellar medium (ISM) in galaxies largely consists of 
hydrogen that occurs in four different phases: the cold neu- 
tral medium (CNM: molecular gas and cold atomic gas), the 
warm neutral medium (WNM: warm atomic gas), the warm 
ionized medium (WIM) and the hot ionized medium (HIM). 
Near the Sun typical temperatures/densities of these phases 
are 80 K740cm~ 3 , 8000 K/0.4cirr 3 , 8000 K/0.2c m - 3 and 
10 6 K/0.003cm- 3 , respectively (Kulkarni & Heiles [19881 
Ferriere 1998 1, which shows that they have about equal ther- 
mal pressures. Physical processes like ionization and cooling 
depend on the structure of the ISM (i.e. diffuse, or cloudy, 
clumped), therefore the volume filling factors of the different 
phases as well as of the gas within these phases are important 
quantities. 

Observations of external galaxies and the Milky Way 
showed that the warm and hot phases are extended and seem 
to occupy most of the interstellar space, whereas the CNM 



fills only a small fraction (Kulkarni & Heiles Q988; Dickey 
119931 Ferriere Tl998> . The WIM consists of dense classical 
H II regions and diffuse ionized gas (DIG) around and in be- 
tween these regions. Also ionized surfaces of cool clouds 
may contribute to the DIG (Miller & Cox [1993 1. Walterbos 
& Braun ( 1994 1 described the wide-spread, diffuse Ha emis- 
sion visible in the spiral arms of M 3 1 as tori around the arms 
with negligible Ha emission in between the arms. 



Reynolds (1991bi obtained an exponential scale height 
of the DIG in the solar neighbourhood of about 900 pc from 
the increase of pulsar dispersion measures perpendicular to 
the Galactic plane with height above the plane. Scale heights 
near 1 kpc were also found by other authors using different 
pulsar samples (Bhattacharya & Verbunt [19911 Nordgren et 
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al.EHU Gomez et al. l2U0H Cordes & Lazio l2~07Bt . The first 
estimates of the volume filling factor of this ionized layer 
were made by Reynolds ( H977L who derived a lower limit 
to the filling factor of ~ 0.1 from the emission measures and 
dispersion measures in the direction of 24 pulsars at Galactic 
latitudes |6| > 5° and a model of the ionized gas layer. In a 
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later paper (Reynolds 1991a) he used the dispersion measures 
and emission measures towards 4 pulsars in globular clusters 
about 3 kpc away from the midplane of the disk. He found 
a mean filling factor of <; 0.2 which represents an average 
value through the full layer. The ionized gas appeared to be 
concentrated in extended clouds of mean density 0.08 cm -3 
occupying about 200 pc along a line of sight towards the 
Galactic pole. 

Maps of neutral and ionized gas of the Milky Way and 
other edge-on galaxies show that bright, dense clouds are 
concentrated to a thin disk that is surrounded by extended, 
diffuse emission from less dense gas reaching large dis- 
tances from the galaxy plane. This indicates that there is 
an anticorrelation between gas density and area filling fac- 
tor with the latter increasing away from the plane. In the 
Milky Way the volume density of the atomic gas decreases 
with increasing distance from the Galactic plane (Dickey & 
Lockman 1990 1, thus the ionized gas density may also de- 
crease and its volume filling factor may increase. Kulkarni 
& Heiles (1988 1 estimated the variation of the local vol- 
ume filling factor, /, and the local mean density within ion- 
ized clouds, n c , with distance from the plane, z, assum- 
ing exponential functions for the local electron densities n e 
and n 2 . They obtained f(z) ~ 0.11 exp(|z|/640 pc) and 
n c (z) — 0.27 exp(— |z|/360pc) cm -3 , scaled to a Galactic 
centre distance of 8.5 kpc. The local volume filling factor 
then increases from 0.11 at z — Opc to 0.52 at \z\ — 1 kpc. 
At | z\ = 400 pc, about half the scale height of the DIG, 
/ ~ 0.21 and n c ~ 0.08 cm -3 , in agreement with Reynold's 
(I1991al results. 

The relationship between the average volume filling fac- 
tor / v and mean electron density within clouds n c , both av- 
eraged along the line of sight, carries important information 
on the structure of the DIG. For example, Elmegreen ( 119981 
1999) showed that turbulence causing hierarchical, fractal 
structure leads to an inverse relationship between f v and n c . 



n c in the DIG. The distances of the pulsars we took from 
the new model of the electron distribution in the Galaxy pre- 



Pynzar dl993> was the first to investigate this relationship, 
not only for the DIG but also for the classical H II regions. 
He used dispersion measures and emission measures towards 
pulsars and observations of recombination lines from dense 
Hll regions. These data yielded f v (n c ) oc n~ . Estimates 
for a mixture of spiral arm and interarm regions by Heiles 



et al. (1996) and Berkhuijsen (1998) are in agreement with 
Pynzar's result (see compilation by Berkhuijsen 1998 ). An in- 
verse relationship between f v and n c is indicated by the near 
constancy of the average electron density along the line of 
sight (n c ) derived from dispersion measures of pulsars with 
known distances. The constancy of (n c ) was first pointed out 
by Weisberg et al. ( 1980); although it is widely used to de- 
rive pulsar distances, its physical significance did not attract 
attention. 

Recently, a new masterlist of pulsar data has become 
available (Hobbs & Manchester 2003 ). We used the disper- 
sion measures of pulsars at Galactic latitudes above |5°| in 
this list together with emission measures from the Wisconsin 
Ha Mapper (WHAM, Haffner et al. 120031 . corrected for ab- 
sorption and for contributions from beyond the pulsars, for an 
improved determination of the relationship between f v and 



sented by Cordes & Lazio (2002). In this sample of several 
hundred pulsars we also searched for variations of f v and n c 
with z in the DIG, determined their exponential scale heights 
and derived the scale heights of the local functions f(z) and 
n c (z). 

The paper is organized as follows: In Sect. 2 we describe 
the derivation of f v and n c and the data used. In Sect. 3 
we explain the statistical procedure and describe the results, 
which are discussed in Sect. 4. Section 5 gives a summary of 
our conclusions. Appendix A contains a glossary of the vari- 
ables used. A table listing relevant parameters of the pulsars 
in the final sample is available on request 1 . 

2. Basic relations and data used 
2.1. Basic relations 

The expressions for dispersion measure, DM, and emission 
measure, EM, towards a pulsar at a distance D (in pc) can be 
written in several ways: 

n e (l)dl = (n c )D = n c f D D = n c L c , (1) 



DM 



cm J pc 

EM 
cm -6 pc 



n 2 Jl)dl = (n 2 c )D = n 2 J D D = n 2 c L c , (2) 



where n c (l) (in cm -3 ) is the electron density at a point 
I along the line of sight, L c (in pc) the total path length 
through the regions containing free electrons (clouds) and n c 
(in cm" 3 ) the average density in these regions which is the 
mean density of a cloud if constant for all clouds along the 
line of sight (see Fig.Q. Furthermore, (n e ) and (rig) are av- 
erages along D and / D = L c /D is the fraction of the line 
of sight occupied by electrons. Note that all quantities with 
overbars are averages along a line of sight. 2 
Combining Eqs. (1) and (2) we find 



/d = 



D 



— 2 



Other useful expressions are: 

_ DM _ EM _ DM 2 — DM 2 
Uc ~ ~L~7 ~ DM ' ~~ "EM" a _ EM~D 



(3) 



(4) 



It is not a priori clear whether the line-of-sight filling fac- 
tor / D is the same as the volume filling factor f v that we want 
to measure. The volume filling factor is the fraction of the ob- 
served volume occupied by the clouds. Generally, / v = c/ D , 
where c is a factor depending on distance, size and shape of 
the clouds. Looking towards a pulsar, the observed volume is 
the volume of the radio beam out to the distance of the pulsar. 
The volume of this cone is Vbcam = jD(D tan |) 2 , where 
9 is the half-power beamwidth in degrees. As 6 is small, 

1 eberkhuijsen@mpifr-bonn.mpg.de 
We keep the conventional forms (n c ) and (n c ), but use overbars 
for other averaged variables to clearly distinguish similar quantities 
and to simplify the notation. 



2.3 Corrections to emission measures 
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Fig. 1. Three ionized clouds along the line of sight towards a 
pulsar P at distance D. The path lengths through the clouds 
- l\, I2 and l 3 - are indicated by thick lines. The total path 
length filled with electrons is L c = l-y + 1% + I3 and the 
filling fraction / D = L c /D. The local electron densities 
within the clouds - n c x, n C 2 and n c3 - give a mean den- 
sity in clouds n c = {n c ih + n C 2h + n c3 l 3 )/L c . The av- 
erage electron density along the line of sight is (n c ) = 
(n c ih + n c2 l 2 + n c3 l 3 )/D. Thus (n c ) =_/ D n c . If / D rep- 
resents the average volume filling factor / v (see Sect. 2.1), 
(«o) = f v n c . 

tan f = 2IJ anc l Vbcam = j^Dd 2 , where d is the linear size 
of the beam at distance D. A thin cloud at distance D and of 
thickness I « D along the line of sight, more extended than 
the beamwidth, occupies V c ~ j^ld 2 , so in this case indeed 
f v ~ / D = Z/d and c ~ 1. Similarly, c ~ 2 for a spherical 
cloud of diameter d equal to the beamwidth, located near the 
pulsar, and c ^ 1 for spherical clouds of diameter I SS e? 1 / 3 . 
As the DIG mainly consists of filaments and sheets (Reynolds 
1991bi, we could assume c ~ 1 and / D = / v . However, an 
extended sheet located close to the pulsar will occupy a larger 
fraction of the beam volume than a similar cloud near the ob- 
server. On average, clouds will be situated about half-way to 
the pulsar, near D/2. Then for a cloud thickness I = fj^D 
the cloud volume is V c = <^Dd 2 [(l + / D ) 3 - (1 - 7d) 3 ] 
andjV = V C /V heam = |](1+7d) 3 -(1-7d) 3 ]. This gives 
c = / v //d = 0-75 for / D ^ 0.15 with a slow increase to 
c = 1 for larger values of / D . 

In most directions there will be several clouds or sheets 
along the line of sight with a total thickness L c that is centred 
near D/2. Therefore we conclude that generally / D ~ f v 
and we will use / v in the remaining of the paper. 

2.2. The data 

In order to obtain realistic values for n c and f v from Eqs. 
(1) to (4) DM and EM should arise in the same ionized re- 
gions in front of the pulsar. For nearby pulsars at low Galac- 
tic latitudes a considerable amount of the observed emission 
measure EM may originate beyond the distance of the pulsar, 
whereas for pulsars at large distances absorption may reduce 
EM. Near the Galactic plane classical H II regions on the line 
of sight towards the pulsar may also influence the results (Mi- 
tra et al. l2003> . To minimize these problems we used pulsars 
at Galactic latitudes \b\ > 5°. 

The dispersion measures of the pulsars were taken from 
the catalogue of Hobbs & Manchester (2003 1. Among them 
are 19 pulsars with accurate parallactic distances (Brisken et 
al. l2002> . We checked their positions for bright H II regions, 
which were found towards 6 of them. We used the 13 distance 



calibrators left for comparison with a much larger sample. Of 
these 13 pulsars 8 are at |6| > 5°. 

The emission measures were taken from the WHAM 
Northern Sky Survey (Haffner et al. l2003> . This data set rep- 
resents the total Ha intensity, integrated between LSR veloc- 
ities —80 and +80 km/s, within a one-degree diameter beam 
centred at a specified position. The intensities and errors 
are given in Rayleighs (10 6 /47r photons cm~ 2 s _1 sr _1 ); 
the sensitivity of the survey is 0.1R. At the Ha wavelength 
1R = 2.41 10 -7 ergcm~ 2 s _1 sr _1 which is equivalent to 
EM = 2.25 cm" 6 pc for T = 8000 K, the mean electron 
temperature in the DIG. A possible increase of T c with in- 
creasing \z\ (Haffner et al. 4999 1 will hardly influence emis- 
sion measures observed through the disk at \b\ > 5°. For the 
interpolation of these data to a pulsar position we used a mod- 
ified quadratic Shepard method for interpolation of scattered 
data in a plane and quadratic error propagation. With the con- 
straint \b\ > 5° the overlap between the pulsar catalogue and 
WHAM survey yields a sample of 320 pulsars. 

In Fig. |2] we plot EM and DM as a function of Galac- 
tic latitude. The sharp rise in both EM and DM at \b\ < 5° 
is due to the classical Hll regions and increasing lines of 
sight. In most cases EM < DM, which indicates that gen- 
erally 7i c < 1 cm~ 3 . For all further analyses we excluded 30 
more pulsars from our \b\ > 5° sample based on the fact that 
they deviate significantly from the outer envelope of the EM 
distribution. We also excluded pulsars located in the Magel- 
lanic clouds and in globular clusters. Our sample at |6| > 5° 
then consisted of 285 pulsars including 8 of the 13 distance 
calibrators. The relevant parameters of the 13 calibrators are 
listed in Tabled 

We determined the distances towards the pulsars in our 
sample from the new model of the electron distribution in the 
Galaxy of Cordes & Lazio ( 2002 ) which is the most realistic 
model presently available. Although its basic structure (i.e. 
thick disk, thin disk, spiral arms) is nearly identical to that of 
the older models of Taylor & Cordes ( 1993 1 and Gomez et al. 
(2001 ), it contains more detail resulting from the larger data 
sets of various kinds on which it is based (i.e. not only on 
dispersion measures, but also on scattering measures as well 
as optical and X-ray observations). Therefore this model per- 
forms considerably better than the earlier models (see Table 1 
in Cordes & Lazio ( 2003 1. It is well defined for the pulsars in 
our final sample, which are at |6| > 5° and within about 3 kpc 
from the Sun in mainly interarm regions (see Fig.[4jl. The sta- 
tistical error in the model distances is about 20% (Cordes & 
Lazio 2002 their Fig. 12), which is much smaller than the 
intrinsic spread in the observed dispersion measures (s.d. of 
factor 2) and emission measures (s.d. of factor 4, see Fig.0. 
This means that the errors in the model distances are negli- 
gible in our statistical study of the relationship between the 
electron density in clouds and filling factors in the DIG, and 
of their variation with height above the Galactic plane. 

2.3. Corrections to emission measures 

The emission measures observed in the directions of the pul- 
sars are not the ideal values for obtaining filling factors and 
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2.3 Corrections to emission measures 



Table 1. Parameters of 13 pulsars with parallactic distances 
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Jname 
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[ cm -3 pc] 


[kpc] 


[kpc] 


[cm -6 pc] 


[ cm -6 pc] 


[cm" 3 ] 




J0332+5434 


145.0 


-1.2 


26.78 


1.03 

—0.12 


0.022±0.003 


10.0±2.9 


46±18 


1.7 ±0.7 


0.015±0.006 


J08 14+7429 


140.0 


31.6 


5.75 


0.433±0.008 


0.227±0.004 


1.79±0.03 


1.2±0.1 


0.20±0.02 


0.065±0.005 


J0953+0755 


228.9 


43.7 


2.97 


0.262±0.005 


0.181±0.003 


3.6±0.9 


1.9±0.5 


0.65±0.17 


0.017±0.005 


J1136+1551 


241.9 


69.2 


4.85 


0.35 ±0.02 


0.327±0.019 


1.5±0.4 


1.1±0.3 


0.23±0.06 


0.059±0.017 


J1932+1059 


47.4 


-3.9 


3.18 


0.331±0.010 


0.023±0.001 


10.5±1.1 


2.9±0.5 


0.92±0.16 


0.010±0.002 


J20 18+2839 


68.1 


-4.0 


14.18 


0.95 ±0.09 


0.066±0.006 


13.2±0.2 


9.5±1.4 


0.67±0.10 


0.022±0.004 


J2022+2854 


68.9 


-4.7 


24.62 


2 3 

z - ;> -0.6 


in +0 08 

U ' ly -0.05 


18.2±1.2 


25.4±3.1 


1.03±0.13 


010+ 005 
u.uiu_ 003 


J2022+5154 


87.9 


8.4 


22.58 


i q +0-3 
Ly -0.2 


0-28 t°;£ 


10.5±2.1 


11.9±2.5 


0.53±0.11 


n 023 +0.006 


J0922+0638 


225.4 


36.4 


27.31 


1-15 Z° fe 


0-68 t^ol 


6.4±1.1 


6.7±1.2 


0.24±0.04 


q 7 +0.025 
'-0.022 


J1537+1155 


19.9 


48.3 


11.62 


1.08 ±g:i« 


0.81 t°o\l 


2.7±0.2 


2.9±0.2 


0.25±0.02 


044+ 007 

U.U++_ .006 


J 1857+0943 


42.3 


3.1 


13.31 


0.79 ±8;?? 


043+ 0015 

U.U+J_ 009 


18.6±2.9 


13.0^;? 


QR+ - 37 


017+ 009 

u - ul ' -0.006 


J1713+0747 


28.8 


25.2 


15.99 


90 + 004 

-0.02 


n qoq+0.017 

.009 


4.3±0.2 


3.9±0.3 


0.25±0.02 


0.072±0.006 


J1744-1134 


14.8 


9.2 


3.14 


35 + 03 

\J.JJ -0.02 


056+ 005 
u.ujo_ 003 


12.0±0.3 


t 7+0.5 
•>■ ' -0.4 


i iq+0.16 


0.008±0.001 



x) Taylor et al. Jl993l . errors are <; 1%; 2) Brisken et al. <2002l 

3 ' Emission measure in front of the pulsar, corrected for absorption (see Sect. 2.3). Errors are one standard deviation 
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Fig. 2. Latitude distribution of observed dispersion measures 
and emission measures towards 744 pulsars in the field of the 
WHAM survey (Haffner et al. I2003i . from which the emis- 
sion measures were taken. The dashed lines are at b = ±5°. 
Some pulsars with exceptionally high emission measures are 
indicated. 

electron densities from Eqs. (1) to (4). On the one hand they 
are too small because of absorption of the Ha emission by 
dust along the line of sight, while on the other hand they are 
too large because of contributions from distances beyond the 
pulsars. We have developed approximate corrections for these 
two competing effects. 

First we correct the observed emission measures for ab- 
sorption. Following Haffner et al. ( 1998 1, the corrected emis- 
sion measure is 

EM C = EM • A with A = e 2 - 2E ( B - v ) , 

where E(B — V) is the colour excess along the line of sight. 
We can estimate E(B — V) from the results of Diplas & 



Savage ( 1994 1 who obtained E(B — V) sin \ b\ as a function 
of distance \z\ above the Galactic plane for a large num- 
ber of stars. As they selected stars seen through the dif- 
fuse ISM, avoiding dust clouds, their result is especially 
suited for correcting the emission measures from the DIG. 
They derived an exponential scale height of the dust layer of 
152 ± 7 pc and [E(B - V) /kpc] = 0.257 ± .010 mag kpc -1 
in the Galactic plane. Integration of this function to infi- 
nite \z\ yields E(B - V) sin |6| = (0.0391 ± 0.0024) and 
A — exp((0.086 ± 0.005)/ sin |fe|). As our sample only con- 
tains pulsars at \b\ > 5°, the largest value of A = 2.68 ±0.16 
(see Fig. |5} and we expect that the observed and corrected 
emission measures contain contributions from all distances 
along the line of sight through the electron layer. The data 
of Diplas & Savage ( 1994 1 show considerable spread around 
their fitted line which they attribute to inhomogeneities in 
the dust distribution. These will cause variations in the ab- 
sorption of the Ha emission and extra spread in the observed 
emission measures. 3 

Since we expect that EM C contains contributions from the 
full line of sight through the electron layer, EM C sin |6| for 
a pulsar represents the absorption-free emission measure of 
the full layer perpendicular to the Galactic plane at the dis- 
tance of the pulsar. It will vary among pulsars due to spread 
in the data and true variations in electron density with posi- 
tion. But if the absorption correction is applicable, the values 
of EM C sin |6| for the pulsar sample should be independent 
of \z\. We checked this for our total sample and found no de- 
pendence on \z\ for longitudes I > 60°, but a significant de- 
crease of EM C sin \b\ with \z\ for 0° < £ < 60° (see Fig.|3}. 
This could be due to unusual absorption along these lines of 



3 The function of A derived from the work of Diplas & Savage 
<1994t at \b\ > 10° is consistent with Eq. 1 given by Dickinson et 
al. ( 2003 l which is based on the dust map of Schlegel et al. (19981 . 
At lower latitudes our values are smaller than those of Dickinson et 
al., because we only consider diffuse dust. 
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sight. The area contains an extended dust cloud with high 
E(B — V)/kpc at 300 pc distance from the Sun (Sherwood 
119741 Lucke I1978t and a chain of dense molecular clouds 
known as the Aquila Rift (Dame et al. 12001b extending to 
b > +10°. As more than 70% of the pulsars at 0° < £ < 60° 
are at |6| < 15°, many lines of sight will cut through dust 
clouds. In a plot of EM against longitude the effect of the 
dust clouds is clearly visible as a deep minimum centred on 
£ = 40° which does not occur in the dispersion measures. 
Such disagreement is not seen at £ > 60°. Because the ab- 
sorption correction taken from Diplas & Savage (1994) only 
applies to the diffuse dust, we had to omit 128 pulsars in the 
range 0° < I < 60°, b > —15° from our sample, leav- 
ing 157 pulsars at |6| > 5° for further analysis. Of these, 14 
pulsars are at 0° < £ < 60°, b < -15°. Only 5 distance 
calibrators of Tableware left in the sample. The distribution 
of these 157 pulsars projected on the Galactic plane is shown 
in Fig. |3 Most of the pulsars in our final sample are located 
in interarm regions and within about 3 kpc from the Sun. 




Height |z| kpc 

Fig. 3. Dependence of the emission measures (in cm~ 6 pc) 
perpendicular to the Galactic plane, corrected for absorption, 
on distance to the plane. The full lines are least-squares re- 
gression lines of log(EM c sin on \z\. If the absorption 
corrections are adequate EM c sin|6| should be independent 
of | a; | as is the case for pulsars at 60° < £ < 360° (b), but not 
for pulsars at 0° < £ < 60° (a) where the slope of the line is 
significant (-0.40 ± 0.09). 

Next we determine which part of the corrected emission 
measure originates in front of the pulsar. The absorption- 
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Fig. 4. Distribution of the final sample of 157 pulsars pro- 
jected on the Galactic plane. Full lines are contours of elec- 
tron density from the model of Cordes & Lazio (2002i. Dot- 
ted lines indicate Galactic longitude. The Sun is located at 
8.5 kpc from the Galactic centre. 

corrected emission measure increases with distance from the 
Galactic plane as 

EM c (z)sin|6| = / n 2 c {z)dz = 7^(0) e" |z|/h dz , (5) 
Jo Jo 

where n 2 (0) is the value of n 2 (z) at \z\ = Opc and h is the 
exponential scale height. The fraction of the total emission 
measure towards the pole coming from the layer below the 
|z|-distance of the pulsar, \z p \, is the ratio ofEq. (5) integrated 
to \z p \ and that integrated through the full layer. Hence, in 
front of the pulsar originates 

EM p sin |6| = EM C sin |b| • (1 - e -l z pl/ h ) . (6) 

At this point we do not know the scale height h, but we do 
know that the mean value of EM C sin b\ for all pulsars should 
equal n%(0)h for the sample considered. In Sect. 4.1 we com- 
bine EM C sin |6| with information on DM sin \b\(z) to evalu- 
ate h, which appears to be in the range 250 < h < 500 pc. 
The minimum value of 280 ±30 pc is best constrained. There- 
fore we used h = 280 pc in Eq. (6). A larger value of h leads 
to a smaller value of the correction factor B = 1 — e~' Zp l/ h . 
We show in Sect. 3.5 that our results are not very sensitive to 
the choice of h in the range 250 < h < 500 pc. 

In Fig. [5] we illustrate the effect of the two corrections to 
EM for pulsars at different latitudes and distances. Especially 
for pulsars at low \z\ (i.e. small D, low b) the correction for 
the contribution to EM from beyond the pulsar dominates the 
absorption correction leading to EM p < EM. 

The final sample of 157 pulsars projected on the WHAM 
survey is shown in Fig.|6] A table collecting the properties of 
these pulsars is available on request 4 . 

3. Results 

In this section we compare the observed as well as the cor- 
rected EM with the DM towards the 157 pulsars in our sample 
(Sect. 3.1), derive the exponential scale heights of (n e ), (n 2 ), 
n c and / v (Sect. 3.2), and present the relationship between 
/ v and n c (Sect. 3.3). 

4 eberkhuij sen @ mpif r-bonn . mpg.de 



3 RESULTS 



3. 1 Dispersion measures and emission measures 




Fig. 6. Positions of the 157 pulsars in the final sample (filled squares) projected onto the distribution of the Ha emission 
(WHAM survey, Haffner et al. l2003> . The map is centred at i = 180°, b = 0° and t = 0° is at the right-hand border. 



As EM and DM do not directly depend on each other, a 
symmetrical-statistical treatment is needed to find their rela- 
tionship. Therefore we calculated the ordinary least-squares 
bisector in the log-log plane (Isobe et al. 1990 1 using the 
modified program code distributed by Feigelson. In all other 
cases we derived the ordinary least-squares regression line of 
Y on X in the relevant units. The quality of a correlation can 
be judged from the correlation coefficient r and the value of 
t resulting from the Student i-test, which gives the probabil- 
ity that the correlation is accidental. The correlation is highly 
significant at the 3cr level if t > 3.1 for a number of indepen- 
dent points N > 100. The values of r, the standard deviation 
in r, <r r , and t are calculated from the equations 

^(Xi-XXYi-Y) 



EM P oc DM 



1.47±0.09 



Naxcry 



1-r 2 
N -2 



t = r 



N -2 
1-r 2 



where X and Y are the mean values of X and Y, respectively, 
and ax and ay the standard deviations in X and Y, The mea- 
surement errors in EM and DM are small compared to the 
uncertainty of 20% in D (Cordes & Lazio 2002 1. Since all 
these errors are much smaller than the intrinsic scatter of the 
points (<; factor 2), all points were given equal weight in the 
fitting procedures. The results are listed in Tabled 

3.1. Dispersion measures and emission measures 

If EM and DM occur in the same ionized regions along the 
fine of sight, we expect EM = n c DM, or EM oc DM b with 
b > 1. To test this we plot in Fig.[7t the uncorrected emission 
measures EM as a function of DM and in Fig. 03 the cor- 
rected emission measures EM p against DM for all 157 pul- 
sars at |6| > 5°. The corrections made to EM have clearly re- 
duced the spread in the data. The bisector fit in Fig. 0> yields 



and with r = 0.76 and t = 15 the 
correlation is highly significant. Because of the even distri- 
bution in latitude of the pulsars and the large fraction of pul- 
sars (72%) at D > 1 kpc, the slope of the bisector has not 
changed after the corrections to EM (see Tabled. Note that 
the components parallel to the Galactic plane, EM p cos \b\ 
and DM cos 1 6|, are even better correlated with nearly the 
same exponent. 

Pynzar H9931 obtained EM oc DM 2 from the lower en- 
velope to the data in his sample. This much steeper slope 
could be due to an increasing overestimate of the emission 
measures towards pulsars at large distances. Overestimates of 
EM could arise from the properties of the spiral arm model 
used by him or from contributions to EM from H II regions 
along the line of sight, especially as two thirds of the pulsars 
he used are at latitudes below |5°|. 

The variations of DM sin |6| and EM p sin \ b\ with height 
above the Galactic plane, \z\, depend on the vertical distri- 
bution of the ionized gas. Figure^ shows that DMsin|6| 
continuously increases up to \z\ ~ 0.9 kpc and then levels 
off. This behaviour reflects that of the more than 100 pul- 
sars with independent distance estimates on which the model 



of Cordes & Lazio (2002i is partly based. A power-law fit 
to the data in Fig. [8^ for \z\ < 0.9 kpc yields DM sin \b\ = 
(22.1±1.3) |z|i p ° c 4±0 05 cm" 3 pc and is highly significant (see 
Table|2j. 

Figure [8j> shows that EM p sin|6| also increases with \z\ 
although not as strongly as DM sin The power-law fit for 
|*| < 0.9 kpc gives EM p sin \b\ = (2.78 ± 0.28) M£ p 4 c 6±an 
cm~ 6 pc with high significance. It is not clear whether 
EM p sin|6| increases beyond \z\ > 1 kpc, because of the 
large spread. The value of EM c sin|6| for the 130 pul- 
sars at \z\ < 0.9kpc is 3.18 ± 0.19cm~ 6 pc. Without ab- 
sorption correction the mean value of EM sin |6| is 2.18 ± 



3.2 Scale heights of (n e ), (rig), n c and f v 
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Fig. 5. Examples of the corrections to the observed emis- 
sion measure as a function of Galactic latitude for an as- 
sumed value EM = 1 cnr 6 pc towards pulsars at distances 
D = 250 pc and 1 kpc. Upper dashed line: correction for 
absorption, A = e 086 / sm The standard deviation in A 
decreases from 0.16 at |6| = 5° to 0.05 at \b\ = 10°, 0.02 at 
|6| = 20° and < 0.01 at |6| > 40°. Dot-dashed lines: correc- 
tion for contributions to EM from beyond the pulsar distance, 
5 = 1- e -l 2 pl/ /! with h = 280 pc. For h = 500 pc B 
is 0. 1 to 0.2 lower than the curves shown. The standard de- 
viation in B is < 0.08 for 250 < h < 500 pc. Full lines: 
corrected emission measure EM p = EM ■ A ■ B. The errors 
in EMp are typically ^15%; they are dominated by the stan- 
dard deviations in the observed EM. See Sect. 2.3 for further 
explanations. 
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Fig. 7. Comparison of emission measure and dispersion mea- 
sure for the sample of 157 pulsars, (a) Observed emission 
measures, EM. (b) Emission measures corrected for absorp- 
tion and contributions from beyond the pulsar distance, EM p . 
The full lines indicate the bisector fits given in TableEl Filled 
circles are pulsars with parallactic distances. 
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0.13 cm 6 pc 5 , thus the mean absorption towards the poles 
is a factor 1.5. 

Comparison of Figs. [8^ and[8j> reveals the much larger 
spread in EM p than in DM. Part of this extra spread arises 
because fluctuations in electron density influence DM lin- 
early but EMp quadratically. The clumpy structure of the dust 
(Diplas & Savage 1994 1 leads to variations in absorption and 
increases the spread in the observed emission measures. Inad- 
equate corrections for absorption and for contributions from 
beyond the pulsars also enhance the spread in EM p compared 
to that in DM. 

3.2. Scale heights of (n e ), (nf ), n c and f v 

The variations of n c and f v with height above the Galactic 
plane are important functions for the understanding of phys- 
ical processes in the warm ionized disk. Therefore we in- 
vestigate the dependencies of (n e ), (n^), n c and / v on \z\ 
in this section. Although these variables represent averages 
along the line of sight 6 , we approximate their |z| -distributions 

5 This value, corresponding to 0.97 ± 0.06i?, is consistent with 
the mean Ha intensity of the survey of Haffner et al. 1 2003] their 
Fig. 13) at |b| = 25°, the mean latitude of the 130 pulsars in our 
sample. 

6 Averages along the line of sight are equal to averages along \z\. 
For example, (n c ) = DM/D = DM sin |6|/|z|. 
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Fig. 8. (a) Dependence of the dispersion measure perpenci- 
ular to the Galactic plane, DM sin on distance above the 
plane, \z\. (b) Same for the corrected emission measure per- 
pendicular to the plane, EM p sin \b\. Full lines indicate the 
power-law fits given in Table[2] Filled circles are pulsars with 
parallactic distances. 

by exponentials and estimate the scale heights H. The rela- 
tionship between H and the scale height h of the local elec- 
tron density n^(z) is given by Eq. (9) in Sect. 3.5. Because 
DMsin|6|(;2) flattens near \z\ ~ 1 kpc (see Fig. [8^), we 
determined the scale heights for pulsars at |z p | < 0.9 kpc. 
There are 130 pulsars in this sample. 

Figure [9] shows the distributions of the various densities 
and / v in the \nY — \z\ plane. The exponential fits are listed 
in Table|2] 

As the model of Cordes & Lazio ( 2002 ) assumes that (n e ) 
is constant, (n c ) is independent of \z\ (Fig.|9^). The spread 
in (n c ) is within a factor of 2 to 3 around the mean; the 
mean value at the midplane (n )o of 0.0205 ± 0.0014 cm -3 
is in good agreement with earlier determinations for the so- 
lar neighbourhood (Reynolds Il991bl Weisberg et al. 1980 1. 
The spread decreases steadily with increasing \z\. This is a 
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3.3 Volume filling factors and n c 



Table 2. Statistical relationships for h — 0.28 kpc 
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1) Bisector in log Y — log X plane. 2) Regression line of log Y on log X. 3) Regression line of In Y on \z\ with slope 
1/6, where b is the scale height H in kpc; +/— means increases/decreases with \z\. All fits are ordinary least-squares 
fits with errors of one standard deviation. Fits with t > 3.1 are significant at the 3cr-level. 4) Error is undefined. 
Units: DM in cm -3 pc, EM and EM P in cm -6 pc, (n e ) and n c in cm -3 , (nl) in cm' 6 , \z\ in kpc, L c in pc. 
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Fig. 9. Electron densities averaged along the line of sight 
and average volume filling factor as function of height 
above the Galactic plane, \z\. (a) average density (n e ) = 
DM sin |6|/|z|. (b) average of the square of the den- 
sity (n 2 ) = EM p sin |6|/|z|. (c) mean density in clouds 
n c = EMp/DM. (d) average volume filling factor f v = 
DM 2 /(EM P D). Full lines indicate the exponential fits for 
\z\ < 0.9 kpc given in Tabled Filled circles are pulsars with 
parallactic distances. 



real 1 2 1 -effect and not just caused by a longer line of sight, 
because in a plot of (n c ) — Dcos \b\ (not shown) the spread 
in (n c ) remains the same up to a distance of 3 kpc, beyond 
which only 13 pulsars are located. 

Figure |9p shows a decrease of (n 2 ) with \z\ correspond- 
ing to a scale height H of about 750 pc. The bottom pan- 
els of Fig. |9] show the variations of n c (left) and f v (right) 
with \z\. The exponential fits show that n c ( | z | ) decreases with 



H ~ 710 pc, whereas / v (|z|) increases with H ~ 670 pc. 
Their midplane values are n c ,o = 0.41 ± 0.06 cm -3 and 
7 v ,o = 0.050 ± 0.009, respectively (see Table EJi. Note 
that the near equality of the scale heights of (n 2 ) and n c , 
and of f v with opposite sign, follows from the relations 
n c = (nl)/(n c ), f v = (n e ) 2 / (n%) and the constant value 
of (n e ) with \z\ (see Eqs. (1) to (4)). The increase of / v with 
\z\ could cause the decreasing spread in (n c ) with \z\ as long 
lines of sight through low-density regions at high | z \ will give 
less variation in (n c ) than short lines of sight through high- 
density regions at low \z\. 

We conclude that up to \z\ = 0.9 kpc (n c ) the scale 
heights H of (n 2 ), n c and / are about 0.7 kpc for this sam- 
ple of pulsars at 60° < £ < 360° and \b\ > 5°. 



3.3. Volume filling factors and n c 

In Fig. El| we present the dependence of the volume fill- 
ing factor f v on the mean electron density in clouds n c 
for all pulsars in the sample (N = 157). The correlation is 
very good: the power-law fit yields f v (n c ) = (0.0184 ± 
0.0011)n c ~ 107±0 03 with a correlation coefficient r = 0.94± 
0.03 and t — 35 (see Table|2jl. It nearly covers two orders of 
magnitude in n c (0.02 - 2 cm -3 ) and / v (0.9 - 0.009). For 
pulsars at \z\ < 0.9 kpc the resulting fit is nearly the same. 
While an inverse correlation between f v and n c is expected 
from the near constancy of (n c ) = f y n c , the combination 
of DM and EM yields the values of n c and / v for which it 
holds. 

The small spread of about a factor 2 in f v for a cer- 
tain n c is remarkable. This is due to two fortunate facts: 1. 
since n c = EM P /DM and J v = DM 2 /EM P D, D and er- 



3.5 Dependence of various relationships on h 
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rors in D cancel in both quantities (see Eqs. (1) to (4)), 2. 
any variation in EM p , which causes the largest spread in the 
various correlations, influences / and n c by the same fac- 
tor but in opposite directions. Therefore points just slide up 
or down along the distribution of points with slope —1. For 
the same reason the corrections to EM have little effect on 
the f v (n c ) relation. Without corrections we find f v (n c ) = 
(0.0193 ± 0.0009)n c - 104±0 03 for N = 157, which holds for 
somewhat higher densities (0.04 < n c < 4 cm" 3 ) than the 
corrected relation. 

The inverse dependence of f v on n c obtained here is 
somewhat steeper than the dependence derived by Pynzar 
iT993l and Berkhuijsen (IT9981 . Their result, / v oc n c ~ ' 7 , 
was based on a large variety of data, including classical H II 
regions, and the pulsar distances were less well known. With- 
out the classical H II regions Pynzar's data indicate a slope 
closer to — 1. We believe that our coherent data set of 157 
pulsars at |6| > 5° more reliably represents the inverse rela- 
tionship between f v and n c for the diffuse ionized gas. For 
mean cloud densities in the range 0.03 < n c < 2cm~ 3 we 
find volume filling factors of 0.8 > f v > 0.01 (see Fig.fTot. 
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Fig. 10. Dependence of the average volume filling factor f v 
on the mean density in clouds n c for the final sample of 157 
pulsars. The full line represents the power-law fit / v (n c ) = 
(0.0184 ± 0.0011)n c - 107±0 03 . The correlation coefficient 
r = 0.94 ± 0.03 and the Student-t = 35 (see Tabled). Filled 
circles are pulsars with parallactic distances. 



3.4. Extent of ionized regions 

As the volume filling factor f v of the DIG increases with 
decreasing mean density n c , also an increase of the extent of 
the ionized regions is expected. This relationship is shown 
in Fig. where we plotted L e as a function of n c for the 
total sample of 157 pulsars. A power-law fit yields L e (n c ) = 
(35 ± 5)tx c ~ - 86±0 06 pc with high significance. The extent 
of the ionized gas along the line of sight towards the pulsars 
ranges from about 35 pc for n c ~ 1 cm~ 3 to about 500 pc 
for n c ~ 0.05 cm -3 . Like in the case of / v (n c ), the slope of 
this relationship is insensitive to errors in EM P . The spread 
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Fig. 11. Dependence of the total path length along the line of 
sight through the ionized regions, L c , on the mean density 
in these regions n c . The full line represents the power-law fit 
for the sample of 157 pulsars (see Table[2j. Filled circles are 
pulsars with parallactic distances. 

in the distribution of L e (n c ) is larger than in that of f v (n c ), 
because L c still contains the distance D. 

Perhaps more interesting than the extent of the ionized 
regions towards the pulsars, which depends on their latitude, 
is their extent in the |z| -direction. This is best obtained from 
the relation L c sin|6|(z) = f v (z)\z\. L e sin|6| grows from 
about 6 pc for \z\ = 100 pc to 53 pc for \z\ = 500 pc and 
about 200 pc near \z \ = 1 kpc. The corresponding mean den- 
sities in the clouds, calculated from n c (z) given in Table 13 
are 0.35 cm -3 , 0.20 cm -3 and 0.12 cm -3 , respectively. 

3.5. Dependence of various relationships on h 

The results presented in Table |2] were obtained using h — 
0.28 kpc to correct the observed emission measures for con- 
tributions from beyond the pulsar distance (see Eq. 6). How- 
ever, in Sect. 4.1 we find that h could be as high as 0.5 kpc. 
To see how h > 0.28 kpc influences the values in Table f2J 
we repeated the fitting procedures for h = 0.36 kpc and 
h = 0.5 kpc. 

The correlations Y(X) corresponding to EM P (DM), 
EM P cos |6|(DMcos |6|) and f v (n c ) remained the same to 
within 1 s.d. from the values in Table [2] while those of 
EM p sin 1 6 1 (DM sin 1 6 1 ) and EM p sin 1 6| ( | z \ ) stayed within 2 
s.d. from the values in Table|2] Only the scale heights H and 
the midplane values Y (h) of (n^), n c and / v changed by 
^ 2 s.d. if h = 0.5 kpc is used. For 0.25 < h < 0.5 kpc 
these values can be found from 

H(h) = 1.8h + 0.24 kpc for (n 2 c ) . (7) 

The scale heights of n c and / v are 0.95H(h) and 0.89H(h), 
respectively. 

Y (h) = 1.04F (0.28) * 0.28/h for (n 2 c ) and n c , 

Y (h) = 1.04Y" (0.28) * /i/0.28 for/ v , (8) 
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4. 1 Vertical structure of the DIG 



where Y" o (0.28) = a in Table|2] 

The general relationship between H and h is given by 



«){z) 



EM p sin|6|(z 
N 



1 /"l z l 

M Jo 



(9) 



As (rig) (z) is not a purely exponential function, the last ex- 
pression is only an approximation. Similar relations hold for 
n c and f v . We checked how well this approximation de- 
scribes the observations by simulating a data set of (n 2 )(z) 
calculated from the integral between \z\ = and 0.9 kpc di- 
vided by \z\ for h = 0.28 kpc, n^(0) = 1cm -6 , and fit- 
ting an exponential to this (noiseless) distribution, exactly as 
we did in Sect. 3.2. We found H = 0.75 ± 0.02 kpc and 
(^e)o = 0.94 ± 0.01 cm" 6 . The errors and the decrease of 
6% in (?ig}o are caused by the curvature of ln(n^)(z) versus 
\z\. Thus the simulation yields the same scale height H as 
derived from the (noisy) observations (see Table |3 and in- 
dicates that the midplane value in Table |2] may be 6% low, 
which is within the observational errors. 



4. Discussion 

In this section we focus on the vertical structure of the elec- 
tron layer. We derive the scale height h of the local electron 
density in Sect. 4.1, compare this with other scale heights in 
Sect. 4.2 and discuss the relationship between f v and n c in 
Sect. 4.3. 

4.1. Vertical structure of the DIG 

The variation of (n c ) with \z\ for all pulsars in our sample 
(N=157) is shown in Fig. [21 Below \z\ ~ 0.9 kpc (n c ) is 
nearly constant (see also Fig.|9^) but beyond this height (n e ) 
starts decreasing. This variation reflects the linear increase of 
DM sin \b\ up to \z\ ~ 0.9 kpc and its levelling off at higher 
\z\ observed for various samples of pulsars at known dis- 
tances (see also Fig.j8ji). Such a behaviour of DM sin \ b\(z) is 
expected if the local electron density n Q (z) decreases expo- 
nentially with a scale height of about 1 kpc (Revnolds ll991bl 
Bhattacharya & Verbunt [T99Tl Nordgren et al. 1 19921 Gomez 
et al. p001 Cordes & Laziop003l which was also adopted 
by Cordes & Lazio ( 2002 ) for their model. However, in an 
analysis of the dependence of DM sin \b\ on \z\ of about 70 
pulsars with independent distance estimates, Cordes & Lazio 



(2003 1 noted that this dependence may be due to a constant 
local electron density up to \z\ ~ 1 kpc and that the electron 
layer may be bounded near this height. Here we give an alter- 
native explanation of the dependence of DM sin |b| on \z\. 

The key observation is the constancy of (n e ) for \z\ 
0.9 kpc. This implies that the local electron density n c (z) is 
also constant and that the local filling factor f(z) and local 
density in clouds n c (z) are inversely correlated. We expect 
n c (z) to decrease with increasing \z\. This is indicated by the 
decrease of emission measure observed for the Perseus arm 
(Haffner et al. 1999 1 and external, edge-on galaxies. Hence, 



f(z) increases with \z\. Beyond \z\ ~ 0.9 kpc DM sin \b\(z) 
flattens and (n c )(z) starts decreasing (see Fig. I12i. We sug- 
gest that this occurs because near \z\ = 0.9 kpc f(z) reaches 
a maximum whereas n c (z) continues to decrease beyond this 
height. DM sin \b\(z) may then be described in the following 
way: 

r\*\ rVA 
DMsin|6|(z) = / n c (z)dz = / f{z)n c (z)dz 



/oe |z|/hf ncoe-'* 



\'^dz 



(10) 



where we assumed that f(z) and n c {z) are exponentials with 
midplane values / and n Ci o, and scale heights hi and h u , 
respectively. 
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Fig. 12. Dependence of the average electron density along 
the line of sight (n c ) = DM sin |6|/|z| on height above the 
Galactic plane \z\ for the sample of 157 pulsars. The full hor- 
izontal line for \z\ < 0.9 kpc is at (n c ) = 0.0213 cm" 3 . The 
continuation to larger \z\ shows the values expected (Eq. (11) 
divided by \z\) for a true scale height h = 280 pc and max- 
imum local filling factor /(zi) = 1 at \z\\ = 0.9 kpc. See 
Sect. 4. 1 for explanations. 



We know that for \z\ < \zi\ ~ 0.9 kpc DM sin \b\(z) = 
(n e )\z\. Hence hf — h n — h, /o«. c ,o = "c(0) = (n c ) and 
f{z) < f{z\) where f{z\) is the maximum value. For \z\ > 
\z\ \ we have 

DM sin \b\(z) = (n c )|zi| + f(zi) / n c , e~^ h dz . (11) 

Thus for very large \z\ DM sin |6| (z) reaches a maximum 

DMsin|6| max = (n c )\z 1 \ + f(z 1 )n cfi he~^ h , (12) 

where |zi| is the height at which the turn-over in 
DM sin \b\(z) starts and f(zi) = f e^ Zl ^ h . 

Similarly, we may describe the emission measure perpen- 
dicular to the Galactic plane as 

EM p sin|6|(z) = / n 2 c (z)dz = / f(z)n 2 c (z)dz . 
Jo Jo 
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Since f(z)n c (z) = f n cfl for \z\ < \z x \ and f(z) = f{zi) 
for \z\ > \z\ , we have 

Az<z x \ 

EM p sin|&|(z) = / n c ,o / n c , e^ z ^ h dz 
Jo 

+ f(zi) [ Z>Z1 nl e-^ h dz. 

Thus the maximum emission measure through the electron 
layer is 

EM p sin|6| max = (n c K, h(l-e-^l/' 1 ) 

+ /(^^e-'W/*. (13) 

This is equal to the absorption-corrected emission measure 
EM C sin |6| discussed in Sect. 2.3. For \z\ > \z±\ we expect 
little increase of the emission measure because the second 
term in Eq. (13) is much smaller than the second term in 
Eq. (12) and even the dispersion measures increase by ^ 25% 
beyond \z\ — 0.9 kpc (see Fig.[8k). Therefore we assume that 
at \z\ — \zi \ the maximum value is nearly reached and reduce 
Eq. (13) to 

EM p sin |6| max = {n c )n cfi h . (14) 

We can now evaluate the scale height h by combining 
Eqs. (12) and (14) with known data. For (n c ) we take the 
mean value of the slightly sloping function (n c )(z), (n c ) = 
/o«c : o = 0.0213 ± 0.0002 cm- 3 . The average value of 
the absorption-corrected emission measures EM c sin|6| = 
3.2 ± 0.2cm~ 6 pc (see Sect. 2.3), which gives n Cy oh — 
150 ± 10 cm -3 pc. From Fig.|8t we obtain DM sin 6| max = 
25 ± 2cm~ 3 pc and z\ ~ 0.9 kpc. We first take f(zi) = 
Jo el 21 1/' 1 = 1 which gives a lower limit to h. For an assumed 
value of h we then calculate n c .o, fo, \zi\, DMsin|6|(zi) 
and DMsin|6| max . We found \zi\ and DMsin|6| max to be 
in the correct range for 250 < h < 310 pc. Repeating the 
procedure for a more realistic value of f(z\) = 0.5 yielded 
a somewhat larger scale height: 320 < h < 400 pc (see 
Table |5Jl. The height |zi| decreases when h increases. The 
lowest value of \zi\ consistent with the data in Fig. [8^ is 
~ 750 pc. Assuming this value gives an upper limit to ft in 
the range 330 < h < 520 pc and a maximum filling fac- 
tor of 0.46 > f{zi) > 0.31. We conclude that the scale 
height of the local electron density in the DIG is in the range 
250 < h Ss 500 pc and that the maximum filling factor 
f(z\) > 0.3. Because the lower limit to h is best constrained, 
we used h — 280 pc to correct the emission measures for 
contributions from beyond the pulsars in Sect. 2.3. How a 
larger value of h changes the results in Table |2] is described 
in Sect. 3.5. 

In Fig.^lthe expected variation of (n e )(z) is shown for 
h = 280 pc and f{zi) — 1. The agreement with the data 
at \z\ > 0.9 kpc is very good. With f(zi) — 0.5 (n c )(z) is 
constant up to \z\\ = 820 pc and then decreases more slowly 
to the same value at \z\ = 2 kpc as for f(zi) = 1. Also for 
\z%\ = 750 pc and h = 500 pc the line is nearly identical to 
those of the two other cases. 

For an electron layer with constant (n c ) up to 1 kpc and 
then a sudden cutoff, as proposed by Cordes & Lazio ( 2003 1, 



(n e )(z) = (n e )(\z\ = lkpc)/|z| for \z\ > 1 kpc. This leads 
to a somewhat faster decrease in (n c ) at high \z\ than is ob- 
served. A more gradual end of the layer seems more realistic 
and might be in better agreement with the data. 

Above we have described that (n e ) (z) will be constant as 
long as the local filling factor and electron density in clouds 
are inversely correlated. This correlation breaks down when 
the filling factor has reached a maximum value of J; 35% 
near \z\ = 0.8 kpc, whereas the electron density in clouds 
continues to decrease to higher \z\. This picture is in good 
agreement with the data in Fig. 1121 and shows that the elec- 
tron layer does not need to be cut off. Why the filling fac- 
tor reaches a maximum at high \z\ is an interesting question. 
Pressure balances must play an important role and model cal- 
culations will be needed to better understand what is going 
on near \z\ = 1 kpc (see also Sect. 4.3.1). 

4.2. Scale heights and column densities 

In the foregoing section we showed that at heights above 
the plane ^ 1 kpc the scale height of the emission mea- 
sure equals the scale height h of the local electron density 
in clouds n c . This happens because the local volume fill- 
ing factor / increases with the same scale height leading to 
n c (z) = (n c ) ~ constant (see Eq. 13). However, at larger 
heights the filling factor remains constant at the maximum 
value f(zi) and the scale height of the emission measure be- 
comes h/2, the scale height of n^. Unfortunately the emis- 
sion measure scale height is difficult to observe and the avail- 
able estimates usually refer to \z\ 1 kpc. 

A lower limit to h comes from radio recombination lines 
observed between I = 332° and 83° near 6 = 0° (Roshi 
& Anantharamaiah 2000 1. These lines originate from low- 
density (1 — 10 cm -3 ) ionized gas in envelopes of normal 
Hll regions. The scale height of this gas is ^ 100 pc. Ha 
observations of M31 (Walterbos & Braun 1994 1 and other 
nearby galaxies (Beckman et al. 2002 1 show that the DIG is 
concentrated around the H II regions in the spiral arms. The 
latter authors argue that electrons leaking from bounded H II 
regions can travel very large distances and could be the origin 
of the DIG at high \z\. 

Reynolds ( 1986 1 estimated the scale height from a com- 
parison of Ha intensities towards the poles and a mean value 
in the Galactic plane. After correction for absorption and 
scattering he obtained a scale height of about 250 pc (cor- 
rected to i?Q = 8.5 kpc), within the range of our value of 
250 < h < 500 pc, although he used a very different method. 

Berkhuijsen et al. ( 1997 1 estimated the scale height of the 
thermal radio continuum emission at 1 .4 GHz, which is unaf- 
fected by absorption, using the data of Reich & Reich ( 1988 1. 
The latitude extent of the thermal emission at galactocentric 
distance R ~ 4 kpc indicates a scale height of about 350 pc 
(after correction for the size of the telescope beam). 

From Ha observations in the direction 125° < I < 152° 
Haffner et al. (0999 » derived a scale height of the emission 
measure of 500 ± 50 pc for distances between 700 pc and 
1750 pc from the Galactic plane in the Perseus arm. Because 
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Table 3. Possible range of true scale height h 
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DMsin|6| max 
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0.536 


0.0397 


904 


19.2 


25.1 




310 


0.484 


0.0440 


970 


20.6 


27.1 


0.5 1 ' 


320 


0.469 


0.0454 


767 


16.3 


23.2 




360 


0.417 


0.0511 


821 


17.5 


25.2 




400 


0.375 


0.0568 


870 


18.5 


27.1 


0.46 


330 


0.455 


0.0469 


750 1 ' 


16.0 


23.1 


0.36 


420 


0.357 


0.0596 






25.0 


0.31 


520 


0.288 


0.0738 






27.1 



assumed value. See Sect. 4. 1 for explanations 



of various line-of-sight effects the observed emission mea- 
sures are dominated by gas at heights smaller than the ge- 
ometrical z-distances and the authors note that the derived 
scale height will be too large. This scale height of 5s 500 pc 
will represent a mixture of the true scale height h expected 
for z < z\ ~ lkpc and h/2 expected for z > Z\, hence 
<, 500 pc < h <> 1000 pc. We note that beyond the Galactic 
solar radius the scale height of the gas increases. Our value of 
250 < h <^ 500 pc for \z\ < \zi \ near the Sun would increase 
to 300 < h <, 600 pc at the radius of ~ 10 kpc of the Perseus 
arm at £ ~ 135° (Merrifield 1 19921 1. Thus if the scale height 
derived for the Perseus arm is more representative for z < z\ 
than for z > z\ it is consistent with our value obtained for 
the solar neighbourhood, otherwise it would be larger than 
expected from the value near the Sun. 

In view of all these data, the value of the true scale height 
derived by us, 250 < h 5s 500 pc for \z\ < \zi\ ~ lkpc 
seems quite plausible for diffuse ionized gas with mean den- 
sities in clouds of 0.03 < n c < 2 cm -3 (see Figs.|9j; and 

In Fig. [HI we plotted the variations of (n c ), n c and / v 
with distance from the plane | z | between and 1 kpc calcu- 
lated from the fitted exponentials in Table|2]which are based 
on h = 280 pc. For h = 500 pc the curves of n c and f v are 
flatter, because n c (0) drops from 0.41 cm -3 to 0.24 cm -3 
and / v (0) increases from 0.050 to 0.084, while the values at 
\z\ = 1 kpc do not change significantly. We also plotted the 
average density along the line of sight of the warm, diffuse 
H I gas, (nni.w) = iV(HI, w)/|z|, derived from the work of 
Diplas & Savage (1994i. This decreases faster than n c be- 
cause the true exponential scale height h is only 195 ± 5 pc 
which corresponds to H = 560 ± 20 pc. Combining this 
with the data for (n e ) yields average degrees of ionization 
of the diffuse gas, I w = (n c )/({n c ) + (nHi, w )), between 
0.053 ± 0.004 at the midplane to 0.24 ± 0.02 at \z\ = 1 kpc. 
It is remarkable that I w (z) and f v (z) are so similar. 

Although in Fig. ^] (nni.w) is everywhere larger than 
(n c ), the local density of warm Hi, nm.w, drops below the 
(constant) local electron density n c — /o«. c o = (n c ) at 
\z\ ~ 560 pc (see Table|4j. This height is close to the result 



of Reynolds ( 1991bi, who found that the total electron den- 
sity becomes larger than the total H I density at \z\ ~ 700 pc. 
This implies that above this height the local degree of ioniza- 
tion I w is larger than 0.5. 

The column densities towards \z\ = lkpc of (n e ) 
and (n H i.w) are N((n c )) = DMsin|6| = (0.68 ± 
0.04)10 20 cm- 2 and iV(HI,w) = (2.19 ± 0.06)10 20 cm- 2 . 
They are similar to the lower values given by Reynolds 
J1991bi . This is understandable because our sample of pul- 
sars traces the diffuse gas at |6| > 5° located in mainly inter- 
arm regions. 



Table 4. Local degree of ionization and volume filling factor 
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0.07 
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400 
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0.17 


0.16 


600 


0.017 
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0.56 


0.34 


0.24 


750 


0.0078 
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0.73 


0.58 


0.32 


800 


0.0061 


3.5 3.2 


0.78 0.76 


0.69 


0.32 


900 


0.0036 


5.9 4.4 


0.86 0.81 


0.99 


0.32 


1000 


0.0022 


6.8 5.9 


0.87 0.85 


1 


0.32 



1 - ) Left-hand values for h = 280 pc, right-hand values for 
h = 500 pc 

n m ,w(z) = 0.366 exp(-|z|/195 pc) (Diplas & Savage ll994l ; 

/w = n c /(n c + n H i, w ); 

n e (|«| < \zi\) = /o^co = 0.021 cm -3 and n c (\z\ > \zi\) = 
f(zi)n c (z), where f(z) = fo exp(|z|//i) and 
n c (z) = n c0 exp(— \z\/h). 

h = 280 pc: |si | = 900 pc, f( Zl ) = I, fo = 0.04, 
n c0 = 0.54 cm" 3 

ft = 500pc: \ Zl \ = 750pc/(zi) = 0.32, / = 0.07, 
n c0 = 0.30 cm -3 
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Fig. 13. Various quantities averaged along \z\ as a function of 
height above the Galactic plane. Straight lines are exponential 
fits given in Table |2] Dots - mean electron density in ionized 
clouds n c = EMp/DM. Circles - average electron density 
(n c ) = DM sin Triangles - average density of warm 

diffuse Hi, (n>Hi,w) = N(HI,w)/\z\ taken from Diplas & 
ravage JT994 ). Plusses - average volume filling factor f v — 
DM 2 /(EM p D). Crosses - average degree of ionization 7 W = 
(n e )/((n e ) + (riHi,w))- EM p was obtained with h = 280 pc. 
With h = 500 pc the curves of n c and f v are flatter than 
shown. See Sect. 4.2 for explanations. 



4.3. Volume filling factors 

In this section we discuss the dependencies of f v on \z\ 
n c derived in Sect. 3.2 and possible implications. 



and 



4.3.1. Dependence of / v on \z\ 

Figures |5Jl and ^] show that the volume filling factor av- 
eraged along the line of sight increases considerably with 
height above the Galactic plane. Looking through the elec- 
tron layer to \z\ = 1 kpc, / v = 0.22 ± 0.04 for h = 280 pc 
and / v = 0.25 ± 0.04 for h — 500 pc, values amazingly 
close to the estimate of f v > 0.2 for the full layer made by 
Reynolds (1991ai based on dispersion measures and emis- 
sion measures towards 4 pulsars in globular clusters. But 
this agreement is accidental. Using the method of Kulka- 
rni & Heiles (19881, Reynolds combined the mean values 
of DM sin |&| and EM sin |6| observed for the 4 pulsars with 
values of (n c )o and (nf)o obtained from other sources and 
did not correct EM for absorption. The combination of data 
from different sources may introduce errors. Furthermore, as 
Reynolds noted, the large intrinsic spread in DM and EM and 
the small number of pulsars make this estimate rather uncer- 
tain. 



The average ionization degree 7 W increases in nearly the 
same way as / v (see Fig.ll3>. If thermal pressure balance ex- 
ists in the DIG both quantities would increase if the temper- 
ature rises with height above the plane. Indications of a tem- 
perature increase with | z | have been reported by Haffner et al. 
J 19991 1 for the Perseus arm at \z\ > 700 pc. Reynolds et al. 



( 1998 1 found that the ionization ratio n e /nm is sensitive to 
temperature. We show this ratio as a function of \z\ in TableE 
Since the filling factor of H I is not known, we calculated the 
ratio of the local mean densities n c (z) = f(z)n c (z) (see 
Sect. 4.1) and nHi,w(z) = 0.366 exp(— |z|/195 pc) (Diplas 
& Savage 1994 1. Between \z\ = and \z\ = 1 kpc the ion- 
ization ratio increases by more than a factor 100, reaching 
1 at \z\ — 560 pc. Thus it seems possible that the tempera- 
ture in the DIG indeed increases with growing distance to the 
Galactic plane. 

In Table 0] we also list the local degree of ionization 
Iw = n c /(n c + riHi,w) an d tne l° ca l filling factors f(z) for 
h = 280 pc and h = 500 pc derived in Sect. 4.1. They all 
increase up to \z\ = 750 pc, where f(z) reaches the maxi- 
mum value of 0.32 if h = 500 pc. At this height I w = 0.73. 
This high value suggests that the maximum value of f(z) 
discussed in Sect. 4.1 occurs when the gas becomes nearly 
fully ionized. The height above the Galactic plane where 
this happens is then indicated by the |z|-distance at which 
DM sin \ b\ (z) starts turning over (see Fig.U). 

4.3.2. Dependence of / v on n c 

The inverse correlation between / v and n c , shown in Fig. 1101 
is the tightest relationship of the 12 cases listed in Table[2] As 
explained in Sect. 3.3, f v (n c ) is independent of distance and 
insensitive to errors in emission measure. 

Because both n c and / v vary with distance from the 
Galactic plane (see Figs. [9]; and|9jl), their correlation could 
just be a |z|-effect. However, an inverse relationship is also 
expected from the constancy of (n e ) with Dcos the dis- 
tance component parallel to the plane: (n e )(0) = 0.0206 ± 
0.0009 cm -3 and the slope of the exponential fit in the 
m(n c ) - Dcos |6| plane of -0.005 ± 0.024, i.e. effectively 
zero. 

Roshi & Anantharamaiah i 



(2001 1 made a survey of radio 
recombination lines of the inner Galaxy at low latitudes. They 
derived filling factors of ^ 0.01 for extended regions of dif- 
fuse ionized gas with densities of 1-10 cm -3 . Their values 
are in good agreement with the extension of our relationship 
in Fig. [10] 

We conclude that the inverse relationship between / v and 
n c seems to hold everywhere in the DIG at least for the den- 
sity range 0.01 <n c < 10 cm~ 3 . This general validity shows 
that it describes a basic property of the DIG, possibly even 
of the entire ISM (Berkhuiisen ll999l l. Cordes et al. dl985> 
estimated a filling factor of io~ 4 - 0±0 - 3 for clumps of about 
1 pc size causing scattering of pulsar signals. The relation- 
ships f v (n c ) and L e (n c ) in Table|2]give f v — 1-3 10~ 4 and 
L c = 0.8 pc for n c = 100 cm" 3 , so these clump proper- 
ties also seem to follow f v oc n c _1 derived for the DIG. 
Furthermore, Gaustad & Van Buren (1993) found a simi- 
lar relationship between filling factor and mean density for 
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clouds of diffuse dust within 400 pc from the Sun: /(rid) = 
(0.056 ± 0.020)7^" 125±0 - 40 for 1 < n d < 10 cm- 3 . For 
dust densities — 1 cm -3 and 10 cm -3 this yields dust fill- 
ing factors of 0.06 and 0.003, respectively, close to the filling 
factors of the DIG for the same densities n c . 

Because errors in emission measure hardly influence the 
/ v (n c ) -relation, we also made some fits without any cor- 
rections to EM. For the sample 60° < £ < 360° (N = 
157) we got the same results as in Table [2] to within the 
lcr-errors, but the range in n c extends to n c = 4 cm -3 , 
twice that of the corrected sample, and / v is correspondingly 
lower. For our original, full sample of 285 pulsars we find 
J v (n c = 1) = 0.0223 ± 0.0010, exponent = -1.02 ± 0.02, 
r = 0.94 ± 0.02 and t = 49. Even for the small sample 
of 13 distance calibrators, 5 of which are at |6| < 5° (see 
Table [TJ, we obtain a highly significant inverse correlation 
in agreement with the larger samples within 3er-errors. Thus 
also without any corrections to EM a relationship between / v 
and n c with the correct slope can be derived for the DIG, but 
the individual values of f v and n c will be wrong by a factor 
depending on the pulsar sample. 

The inverse relationship between / v and n c derived here 
refers to mainly interarm regions with relatively low average 
electron density, (n c ) = 0.021 cm" 3 . In spiral arms, where 
(n c ) seems 2-3 times higher (Cordes & Lazio 2002 1, corre- 
spondingly higher filling factors are expected if the slope of 
the relation is near —1. 

Which physical process or processes could cause the 
tight, inverse relationship between f v and n c ? Without go- 
ing into details we mention two possibilities: 

— Thermal pressure equilibrium in the DIG. If the tempera- 
ture increases, the clouds expand and the electron density 
in clouds decreases. In Sect. 4.3.1 we gave some evidence 
for an increase of the temperature with \z\. 

— Turbulence in the ISM causing fractal structure. Much of 
the ISM appears to have fractal structure which leads to 
7 V oc n c -1 (Fleck [T996l Elmegreen [T998l[T999t . 

Determinations of the f v (n c ) -relation for other phases 
of the ISM will be needed to see how widely it is valid. 
The physical processes causing the inverse relationship / v oc 
n^ 1 could be studied in simulations modelling the ISM. 

5. Summary and conclusions 

We have used dispersion measures of pulsars (Hobbs & 
Manchester I2003i . distances from the model of Cordes & 
Lazio ( 2002 1 and emission measures from the WHAM survey 
(Haffner et al. 2003 1 for a statistical study of several proper- 
ties of the diffuse ionized gas (DIG) in the Milky Way. The 
emission measures were corrected for absorption and con- 
tributions from beyond the pulsar as described in Sect. 2.3. 
To avoid regions with excessive absorption and contributions 
from bright H II regions, we selected pulsars at Galactic lati- 
tude \b\ > 5°. The final sample analyzed contains 157 pulsars 
in the longitude range 60° < I < 360°. Their distribution 
projected on the Galactic plane is shown in Fig.|2] The statis- 
tical results are given in Table[2] 



Our main conclusions are summarized below. 

1 . Dispersion measure DM and corrected emission measure 
EM p are well correlated (Fig.^b), indicating that they probe 
the same ionized regions. This allows us to derive and com- 
pare the average densities along the line of sight (n e ) and 
(rig), the mean electron density in clouds along the line of 
sight n c , the volume filling factor of these clouds f v and the 
total path length through the ionized regions L c (see Eqs. (1)- 
(4)). 

2. The total extent of the ionized regions perpendicular to 
the Galactic plane increases from about 6 pc towards \z\ ~ 
100 pc to about 220 pc towards \z\ = 1 kpc, while the mean 
cloud density n c decreases from about 0.35 cm" 3 to about 
0.10 cm" 3 (Sect. 3.4). 

3. Below \z\ ~ 0.9 kpc (n c ) is essentially constant with 
a maximum spread of a factor 2 about the mean value 
of 0.021 ± 0.001 cm" 3 (Fig. [HJ. As (n e ) = 7 v n c , this 
means that f v and n c are inversely correlated. While n c 
decreases with increasing \z\ (Fig. |9p), f v increases with 
\z\ (Fig. |9ji). We derived the relation / v (n c ) = (0.0184 ± 
0.0011)7I C - 107±0 03 , which is very tight. Figure flOl shows 
that it holds for the ranges 0.03 Ss n c < 2 cm" 3 and 
0.8 i; f v > 0.01. The inverse dependence of f v on n c could 
mean that the DIG is in thermal pressure equilibrium or that 
it has a turbulent, fractal structure. 

4. The linear increase of DMsin|6| with \z\ up to \z\ ~ 
0.9 kpc indicates that the local electron density n c (z) is 
constant, so also locally n c (z) and f(z) are inversely cor- 
related. We suggest that the turn-over of DMsin|6| near 
\z\ = 0.9 kpc occurs when f(z) reaches a maximum value, 
whereas n c (z) continues to decrease towards higher \z \ caus- 
ing a flattening of DM sin |&|(z). 

5. Using this interpretation and the mean of the absorption- 
corrected emission measures through the full layer 
EM c sin|6|, we derived a true scale height of n c (z) 
and f(z) between 250 pc and ~ 500 pc for maximum values 
of f(z) between 1 and 0.3, respectively, near \z\ ~ lkpc. 
Since n^(z) = /(z)n^(z) and n e (z) — f(z)n c (z) = con- 
stant, nl(z) has the same scale height as n c (z) (Sect. 4.2). 
Beyond \ z\ ~ 1 kpc, where f(z) ~ constant, the scale height 
of nl(z) becomes a factor of 2 smaller than that of n c (z). 

6. The variation of f v (z) is similar to that of the average de- 
gree of ionization of the warm, atomic gas, I w (z) (Fig. ll3l >. 
Towards \z\ = lkpc/ v = 0.24±0.05 andl w = 0.24±0.02. 
The increase of I w with \z\ may point to an increase of 
the temperature in the DIG with increasing distance to the 
Galactic plane. The local function I w (z) reaches values near 
unity at \z\ = lkpc (Table 0. If f(z) reaches a maxi- 
mum value when the gas is nearly fully ionized, then the \z\- 
distance where this occurs is indicated by the turn-over point 
of DMsin|6|(z). 

As the tight inverse relation between / v and n c is inde- 
pendent of \z\ or distance along the plane, it seems a basic 
property of the DIG. Simulations modelling the DIG may be 
able to reproduce this relationship. 
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Appendix A: Glossary 



A 


absorption correction applied to EM 


B 


correction for emission measure from beyond the pulsar 


D 


distance to pulsar 


DM 


observed dispersion measure 


EM 


observed emission measure 


EM C 


EM corrected for absorption, EM C = EM • A 


EM P 


EM corrected for absorption and contributions from beyond the pulsar, EM P = EM • AB 


/ 


local volume filling factor 


/d 


fraction of line of sight containing electrons 


/v 


average volume filling factor, average fraction of radio beam up to pulsar containing electrons 


H 


exponential scale height of quantities averaged along the line of sight 


h 


true scale height, exponential scale height of local quantities 




local degree of ionization of warm H I 


/w 


degree of ionization of warm H I averaged along the line of sight 


ie 


total path length through ionized regions along D 


n c 


local electron density in a cloud 


n c 


mean electron density within clouds 


n 


local electron density, n c — fn c and n\ — fn\ 


K) 


average electron density along the line of sight, (n c ) = f v n c 


(™c> 


average along the line of sight of the square of the local electron density, (nl) — f v nl 


«HI,w 


local density of warm H I 


(«hi,w) density of warm H I averaged along the line of sight 


z 


height above the Galactic plane 


z p 


height of the pulsar 


Z\ 


height at which / reaches its maximum value 



